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1. Introduction 

Fractional Brownian motion (fBm) is frequently used to model short- and long- 
range dependence. By definition, fBm with Hurst parameter H € (0, 1) is a 
centered Gaussian process { , t > } with the covariance function 

For H > 1/2 fBm has a property of long-range dependence, for < 1/2 it has 
short-range dependence and, in fact, is counterpersistent, i.e. its increments are 
negatively correlated. For H = 1/2 fBm is a standard Wiener process. 

Two important properties of fBm are stationary increments and self-similarity. 
However, this properties restrict applications of fBm, so recently so-called mul- 
tifractional processes gained huge attention. Multifractionality can consist both 
in dependence of memory depth and regularity of process on the time instance 
and on the time scale. In this paper, we are dealing with the latter kind of mul- 
tifractionality, where the properties of the process depend on the size of time 
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interval, on which the process is considered. In other words, we are consider- 
ing processes, which are not inherently self-similar. A simplest approach is to 
consider a linear combination of fBm's with different Hurst parameters. 

Here we will concentrate on the case where we have only two fBm's and one 
of them has Hurst parameter equal to 1/2, simply put, it is a Wiener process. 
So we consider a process 

Mt" ^aB^ + bWt,t>0 (1) 

where a and b are some non-zero coefficients. Such mixed models and their 
applications where considered in many papers, see [1, 6, 11, 16, 17]. 

The main aim of this paper is statistical identification of model (1), i.e. the 
statistical estimation of the model parameters. The principal attention will be 
given to the estimation of H, though we will give some comments on estimation 
of a and b. 

In a "pure" fBm case, there exist several methods to estimate the Hurst 
parameter, an extensive overview of which is given in [7]. The most popular 
methods are based on quadratic and, more generally, higher power variations 
of the process. A huge literature is devoted to such questions, we will cite only 
few: asymptotic behaviour of power variations and, more generally, of non-linear 
transformations of stationary Gaussian sequences is studied in [4, 10, 13, 20], 
and stochastic estimation for £Bm and multifractional processes with the help 
of power variations, in [3, 8, 9, 12, 14, 21]. Weighted power variations serving 
similar purposes for stochastic differential equations driven by fBm, were studied 
in [18, 19]. 

In this paper we will use power variations of process (1) for estimation of 
parameters. Related questions for mixed models were already considered in [5, 
11, 15, 23], but we use different techniques here. We remark that, in contrast to 
the pure fractional case, there is no self-similarity property in mixed model (1), 
so we cannot directly apply the results of [4, 10, 13, 20] on asymptotic behaviour 
of sums of transformed stationary Gaussian sequences. For this reason we need 
to study asymptotic behaviour as n — )■ cxd of "mixed" power variations of the 
form 

n-l 

E i^ik+D/n - Wk/nY {Bfk+l)/n " ^f/n) ' , (2) 
fc=0 

involving increments of independent fBm and Wiener process W, where 
p > 0, r > are fixed integer parameters. We obtain that the behaviour of such 
power variations depends on parity of p and r and on the value of the Hurst 
parameter H. 

The paper is organized as follows. Section 2 contains results on the asymptotic 
behaviour of mixed power variations. These results are used in Section 3 to 
construct strongly consistent estimators of parameters H,a,b in model (1) and 
asymptotic confidence intervals for H. Finally, in Section 4 we present simulation 
results to illustrate quality of the estimators provided. 
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2. Asymptotic behaviour of mixed power variations 

Let W = {Wt,t > 0} be a standard Wiener process and = {B^,t > 0} 
be an independent of W fractional Brownian motion with Hurst parameter 
H E (0, 1) defined on a complete probability space (il, J^, P). 

For a function X : [0, 1] — )■ M and integer n> 1, fc = 0, l,...,n— Iwe denote 
A^X = — Xf^^n- In this section we will study the asymptotic behavior 

as n — > oo of the following mixed power variations 

n-1 
k=0 

where p > 0, r > are fixed integer numbers. Thanks to self-similarity of 
and W, the sequence {{n^ A^B^ ,n^/'^A^W),0 < k < n — l} is equivalent in 
law to {(Cfc, Cfc), < fc < n — 1}, where {^k, A: > 0} is a sequence of independent 
standard Gaussian random variables, {Cfc,fc > 1} is independent of {£,k,k > 1} 
stationary sequence of standard Gaussian variables with covariance 

PH{m) = E [CoCn] = E - B,^)] = ^ (|m + if" + \m- lf")-\m\' 

Therefore, by ergodic theorem, 

fc=0 

a.s., where for integer m > 

M„ = E [iV(0, 1)™ ] = (m - 1)!!1„ even 

is the mth moment of the standard Gaussian law. So it is natural to study the 
following centered sums 

Sn^'-'' = E i^''"^'^' (^^^)' i^'^B^y - ^'p^'^) ■ 

k=0 

The following theorem summarizes limit behaviour of S^ ^'^. We remark 
that some (but not all) of results can be obtained from the limit theorems for 
stationary Gaussian sequences of vectors, see e.g. [2]. However, we believe that 
our approach (using one-dimensional limit theorems) is more accessible and 
leads quicker to the desired results. 

Theorem 2.1. If p and r are even, then 

• forHe (0,3/4) 
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where 

2 V / / \2l 2 / 2\ 



/or if = 3/4 



3/4,p,r 



-^(0' <^3/4,r/"p + <^r,p). " ^ OO, (5) 



V n log n 

where (73/4 = 3r() — l)/4; 

• for He (3/4, 1) 

where C,H,p,r is a special "Rosenblatt" random variable. 
If p is odd, then for any H G (0, 1) 

//p is even and r is odd, then 

• forH€ (0, 1/2] 

w/iere 

(r-l)/2 2 

""'^^ 2. (2/ + l)!((r-2/-l)!!)2 2- ^^^^^ ' 
;=o ^ ' ' m=— 00 

• forHG (1/2, 1) 

n-«5^'f ^ 7V(0, n ^ 00. (9) 

Proof. We study different cases in the same order as they appear in the formu- 
lation. 

Assume first that p and r are even. The principal idea in this case is to rewrite 
mixed power variation as 

SH,p,r Q' _]_ Q" 

where 

n-l 

fe = .^Q^ 

n— 1 ^ ^ 

fe=0 

and to apply known results concerning asymptotic behaviour of S[[, since it 
contains only fractional Brownian motion, and consider S'^ conditionally on the 
fractional Brownian motion. In what follows we realize this idea. 
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For e (0,3/4) write 
According to [4], for r even, H € (0,3/4), 



n— 1 

- V2 {^IB^y ^ Nio, (h) 



k=0 



as n — >■ cx). Consequently, n ^/^S*" A^(0, crj^^^p), n ^ oo. Further, 
where 

_! n — 1 



k=0 
n-1 



1/2 



(12) 



fc=0 



Since is uniformly continuous a.s., niaxi<fc<„ A/j_„ — )■ 0,n — > (X) a.s. Thus, 
taking into account independence of and W and an evident fact that 
vari?„ = 1, we get by CLT that the conditional distribution of Rn given 
converges to standard normal distribution as n — c» a.s. Further, from the 
ergodic theorem 



a.s., hence 



'"-^ Y^iKB^ f ^ ^Jt2r, n^oo 

k=0 



n oo a.s. 



So by Slutsky's theorem, the conditional distribution of S*^ given B^ converges 
to iV(0, ap j.) a.s., that is, for any t G R we have 



J e'*^P(S'; e dx I B") 



(13) 



a.s. as n — > oo. Now write 



Jt{S'„+S'^ 



B 



H 



H\ JtS'' 



e"^P(S'; G dx\B") e 



whence 



< £/! + £'2 
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where 



El = 



f e'*"P(5; G dx\B")-e- 
Jr 



0, rt — > oo 



by (13) and dominated convergence; 



Jts'' 



+2 2 2 /o 



— ?> 0, n — >■ oo 



by (11). It follows that 

as required in this case. 

In the case where if = 3/4, we have by [4] 



n— 1 



fc=0 



as n — >■ cx), whence (5) can be deduced using the same reasoning as above. 
For H £ (3/4, 1), write 

where i?„, fc„, 5*^' are defined above. As before, i?„ N{{), 1) conditionally given 
as rt — >■ oo a.s. However, this time 7i^~^-^+''^fc„ ^> 0, n — >■ oo a.s., since 
has a finite limit and rt'l'^-'^i^ — > 0, ?i — > oo. Therefore, n}~'^^ S'^ 
0, n — > oo. Further, according to [10], see also [13, 20], 



^-2H Qii 



Sn CH,p,r, ri OO, 

where Cff,p,r is a "Rosenblatt" random variable. Thus, we get (6) using Slutsky's 
theorem. This finishes the case where p and r are even. 

Now assume that p or r is odd. In this case S^'^'^ has a form 



S^^P.- = J2 {AlWf {^IB") 



k=0 



Write 



where i?„ and fc„ are defined above, Z„ = '^'^^^^^Mp X]fc=o i^k^^Y ■ 
fore, given i?^, n'"^^^/^i?„fc„ A^(0,(Tp ,,), n — ^ oo, a.s. 

Now if p is odd, we have ~ irrespective of value of H, whence (7) 
immediately follows. 

Further, assume that p is even and r odd. For ii e (0, 1/2], we have by [4] 

n-l 

J2 (A^B^y ^ iV(0, aj,A n ^ oo. 

k=Q 
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Therefore, Z„ N{0, fi^cr^ n -> oo. Arguing as in deriving of (4), we get 
(8). For H e (1/2, 1), it folfows from [10] that 

n-1 

n('-i)^ ^ {AlB^Y ^ N{0, ^jtl^,), n ^ ex., (14) 

fc=0 

so ni/2-^Z„ ^ N{Q, AipAir+i)> oo But 

whence (9) follows. The proof is now complete. □ 
3. Statistical estimation in mixed model 

Now we turn to the question of parametric estimation in the mixed model 

Aff = aSf + bWt, i € [0, T] , (15) 

where a, h are non-zero numbers, which we assume to be positive, without loss 
of generality. Out primary goal is to construct a strongly consistent estimator 
for the Hurst parameter if, given a single observation of . 

It is well-known (see [6]) that for H e (3/4, 1) the measure induced by in 
C[0,T] is equivalent to that of hW . Therefore, the property of almost sure con- 
vergence in this case is independent of H. Consequently, no strongly consistent 
estimator for H based on a single observation of exists. 

In this section we denote A^X — XTk/n ^ ^T(fe-i)/n ^nd 

n-1 
k=0 

Also for a real sequence {a„,n > 1}, non-zero sequence {bmn > 1}, and a 
sequence {^„, n > 1} of random variables, we will abbreviate 

6, « a„ + bnN{0, 1) 

whenever the following two assertions hold: 

• either ^„ ~ a„ as n ^ cx) a.s. or a„ = 0; 

• (^n — a.n)/bn =^ N{0, 1) as n — oo. 

3.1. Statistical estimation based on quadratic variation 
Consider quadratic variation of , i.e. 

n-1 

Kf^^ = 5: (A^M^)^ 

k=0 

The following lemma contains information on asymptotic behaviour of V^^'^, 
which is needed to construct estimators. 
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Lemma 3.1. 1. For H e (0, 1/2), 

V^-^ « a2y2ff^i-2if ^ ^2 J. ^ a^T-'"n^'''-^"<JH;2N{Q, 1). (16) 

2. ForH£ (1/2,1) 

'2 ^ ^ ^2j,2H^i-2H ^ y262Tn-i/2Ar(o, 1). (17) 

3. For H ^ 1/2 

V;^'2 p. (^2 ^ ^2~)j. ^ y2(a2 + b^)Tn-^'^N{0, 1). 

Proof. Write 

^2 ^ „2^H,o,2 2a6V„^'i'i + fo^V^f >2'0. (18) 
By law of large numbers 

yH, 2,0 n^oo. 

Therefore, we have 

yH,2,o ^j.^ V2Tn-^/^N{0, 1) 

by the classical CLT. 

This gives immediately the result in the case H = 1/2, since in this case 

= Va^ + b^Wt, where W is a standard Wiener process. 
Further, by ergodic theorem and (11), for H G (0,3/4) 

'0^2 ^ nl-2Hj.2H ^ n'/2-2"T^"aH,2NiO, 1), 

for H = 3/4 

yH,0,2 ^ ^3/2^-1/2 ^ 3^.3/2,^-1^1^^(0^ 



and for H e (3/4, 1) 



where (^h is a Rosenblatt random variable. 
Finally, from (7), 

yH,l,l^„-ffrff+l/2^(0,l). 

The required statement now follows by comparing the orders of terms in 
(18). □ 

As we see, the asymptotic behaviour of V^'^ depends on whether H < 1/2 or 
not. Let us consider these cases individually (we omit the case H — 1/2, which 
is of minor interest). 
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3.1.1. He (0,1/2) 

For H e (0,1/2) we have 

a.s., which allows to estimate the Hurst parameter H consistently. Precisely, we 
have the following result. 

Proposition 3.1. For H € (0, 1/2), the following statistics 



and 

=2 ,,Ha 



1 /. v"-^ 



are strongly consistent estimators of the Hurst parameter H . 
Proof. Write from (16) 

l0g2 V^,^ = l0g2 (a2r^^2'=(l-2H)^ + log^ ^ b^rj.l-2H^-k(l-2H) ^ 

where the sequence {Cfc, fc > 1} is bounded in probability. Hence we have 
log2 V^^^ - 2 log2 a + 2H \og^ T + (1 - 2H)k 

(20) 



+0(2~''-(i-2«)) + Op{2-^'^),k 



— ?► 00, 



where O is the standard big-0 notation, while = Op{ak), fc — > 00, means 
that the sequence {^k/ctk, fc > 1} is bounded in probability. In particular, 

log2 V^'^ - 2 log2 a + 2H log^ T + (1 - 2H)k, fc cx), 

whence the result immediately follows. □ 

Remark 3.1. At the first sight, there is no clear advantage of Hk or Hk- But a 
careful analysis shows that Hk is better. Indeed, from (20) it is easy to see that 

fr „ log2 a + log2 T 1 

Hk ~ H , h o(k ),fc— >oo, 

fc 

while 

Hk = H + 0(2-'=(i-2H)) ^ 0^(2-/^/2)^ fc ^ 00. (21) 

Now it is absolutely clear that Hk performs much better (unless one hits the 
jackpot by having aT^ = !)• 
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Asymptotic confidence intervals for H and consistent estimates for 

a and h Asymptotic expansion (16) allows to construct confidence intervals 
for H. In view of (16), we should distinguish two cases: H G (0, 1/4) and H G 
[1/4,1). 

In order to proceed, we should first give a consistent estimate for a. The 
previous remark enables us to do this. 

Proposition 3.2. For H e (0,1/4), the statistic 

is a consistent estimator of . 
Proof. Observe that 

in probability, since k{Hk — H) — > 0, fc — > oo, in probability by (21). Therefore, 
we have the desired statement. □ 

The construction of asymptotic confidence intervals is different in cases H € 
(0, 1/4) and H € [1/4, 1). Let us first concentrate on the case H e (0, 1/4). In 
this case 

^ ^2rj.2H^i-2H ^ a^T^"v}'^-^"aHaN{Q, 1). 
Now in order to construct a confidence interval for H , observe that 

where the constant cr^ can be computed explicitly in terms of H , however, the 
expression is extremely cumbersome; thus we give only a less explicit expression 
in terms of : 

/ oo \ 1/2 



fe=l 



where 

PH.k - E 



2i-^(sn'-(sf/2)"-(sf-sf/2 



2 / „ „ \ 2 



X ( 2^ 2^ {B^_^_^ - B")^ - {b^^^i^ - B"^ - {b"^-^ - 
Now, it follows from Slutsky's theorem that 
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whence the construction of asymptotic confidence interval is obvious. However, 
this leads to solving highly complicated equation involving quite unpleasant 
function cr^. A better idea is to replace H with Hk'- 

2^(25.-1/2)^ ^ ^ 2k/2i ^2" ^2^^^ ^ 1) ^ ^ 

Hence, we have the following result, where Za is the standard normal (1 — a/2)- 
quantile; for well-definiteness we denote 




log2+ X = 

Proposition 3.3. Let H e (0, 1/4). The interval [-ffmin , ^^max] , where 

Hnun = ^ (l + l0g2 Vf,'^ - log2 {vf,!,^, + Z ^Cj' ^ V^,^'' 2' ' ^fj , 
ffmax = i^min V ^ (l + log2 V^,^'' - log^^ {v^,f, - Z ^a' ^ V^,^'' 2' ' ^)) , 

is an asymptotic confidence interval for H of level I — a. 

Remark 3.2. As it is known from [22], for He (0, 1/4) the measure induced by 
in C[0,T] is equivalent to that of aB^ . So the results are essentially the 
same in this case as for fractional Brownian motion alone. Another important 
consequence that in this case it is not possible to estimate b consistently. 

Now let us move to the case H e (1/4, 1/2). (We omit the case H = 1/4 for 
two reasons: considered similarly to H E (1/4, 1/2), and it is hard to distinguish 
it statistically from the case where H ^ 1/4.) The construction of confidence 
interval is more involved in this case, since the second term in (16) cannot be 
omitted. However, the idea is the same: one needs to remove non-random terms 
from the (19). To be more precise, we have the following result, where 

jjH,2 TrH,2 T^//.2 

and the asymptotic variance is 

Coo 
k=l 



CffSo,! - (^C^f + Ij I Sq +Si/2)+So + *1/4 + -^1/2 + *3/4 




( 1 / , 1 \ 1/2 , 1/2 \ , 1/4 
: [cHSk - {CH + 1) [S^ + ) + 



1/4 , 1/4 1/4 \ 

''fc+1/4 ' ''fc-l-1/2 ^k+3/4 J 



With 4 = (i?,^,,-Bf)', CH-21 



-2H 
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Proposition 3.4. Let H E [1/4,1/2). The interval [i?mini ^^max]; where 

i^max = ifmin V (l + log2+ (-f/f - log2+ (-(/f/^ - ¥^^^2'''/^ 



2 

is an asymptotic confidence interval for H of level 1 — a. 
Proof. It follows from [4] that 

r,l-2HTrH,2 /-, , r,l-2H\TrH,2 , TrH,2 ^ 2 rp2ffr,(l/2-2H)fe AT/n 1^ 

2 V^k' -(1 + 2 )V^k+i + ^2''+^ ~ '^hT 2^' ' N{0,1). 
Thus, by Slutsky's theorem, we have the weak convergence 

2^/'- ^-^ ^ y H,2 ^ ^(0. 1)' ^ (23) 

'^k^2'=' 

whence the result follows. □ 

Now we turn to estimation of b. In order to proceed, we will need another 
consistent estimator of H. 

Proposition 3.5. For H e (0,1/2), the statistic 

is a consistent estimator of H , moreover, 

=H + Op{2'^l^), k-^oo. 

Remark 3.3. Despite Hj^ has asymptotically a better rate of approximation 
that Hk for H G (1/4, 1/2), we still do not recommend to use it, as the asymp- 
totic variance is high; it's practically useless for k < 10. 

Proof. From (19) and [4] it follows that 

C/f « a^T^" {2"-^" - 1) 2'=(i-2^^) + a2y2H2fe(i/2-2H)^^'^(o, 1) 

with some constant cr^, which is of no importance now. In particular, U^'^ is 
eventually negative, and we can write 

log2(-C/f - 2 log2 a + 2H log^ T 
+ log2 (1 - 2i-2ff) + _ 2i/) + Op(2-*/2). 



Thus, 



log2 =2H-1 + Op(2-'=/2), fc ^ oo, 

whence the assertion follows. □ 
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Now we are able to construct a consistent estimator for b. However, the 
variance of this estimator is so huge that it is practically useless (see Section 4). 

Proposition 3.6. For H £ (1/4,1/2) the statistic 

" (2i-2i^l'' _ l)T 

is a consistent estimator of . 
Proof. Define 

b- 



It easily follows from (16) that bj, b , k ^ oo. So it is enough to show that 

l-bl 



6| — 6| 0, fc ^ oo. To this end, write 



bl-b' 



T(2) 



(2i-2//r' - i)r 

Obviously, the second term converges to zero. To study the first one, write 
whence the consistency of b1 follows. □ 



3.1.2. H e (1/2,3/4) 

Now we move to the case e (1/2, 1). In view of (17), both Hf^ and converge 
to 1/2 for H e [1/2, 1), so they are not suitable for estimating H . 

In order to estimate 7f, the trick is to cancel out the first term of asymptotic 
expansion (17). The most natural way to accomplish this is to consider 

U^'^ - V^.-^, = 2 l^t;- ^ty.M" . (25) 

By (17), 

[/f .2 ^ a2y2H2(l-2ff)fc (1 „ 2l-2^) , fc ^ OO (26) 

a.s., whence we can easily construct estimators exactly as in the case H G 
(0,1/2). 
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Proposition 3.7. For H e (1/2,3/4), statistics 



and 



are strongly consistent estimators of the Hurst parameter H . 

We will see in Section 4 that performs very poorly, and H^. performs much 
better, despite having a worse asymptotic rate of convergence. For the values of 
H close to 1/2, both estimators are quite poor, and it is better to estimate H 
by regressing several values of log2 U^'"^ on k. 

Asymptotic confidence intervals for H and consistent estimates for a 

and h To construct a confidence interval, it is not hard to check (since we are 
not going to use this result, we will not give a proof) that 



2l-2H^H,2 _ ^H,2 ^ 52y^23-4ff + 1 2-'''^-^'^ N {0, 1). 

In contrast to the case H e (0, 1/2), where one has to replace H by its estimate, 
this result can be used immediately to construct asymptotic confidence intervals 
for H (though one needs to replace b^T by V^''^). Indeed, one can write 

where Za is the standard Gaussian (1 — Q;/2)-quantile. Thus, the construction 
of asymptotic confidence interval is reduced to solving quadratic inequalities 
and taking logarithms. Unfortunately, even for such large values as a = 0.1 
and k — 20 these confidence intervals are useless: approximation (27) is very 
imprecise, and, what is more dramatic, they usually contain the whole interval 
(1/2,3/4).^ 

Estimating b is quite straightforward now: it is easy to verify that 

bl = T-'V^^' 

is a consistent estimator b^. We did not manage to construct a reasonable 
variance-based estimator for a. Similar to estimating b in the case H E (0, 1/2), 
this requires to construct an estimate of H with approximation rate o(fc^^). 
However, none of the estimators we presented has this order. One should take 
a much finer partition for this, moreover, the number of points should be de- 
pending on H, which is practically impossible. 
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3.1.3. H e (3/4,1) 

As we have already mentioned in the beginning of this section, in this case it is 
impossible to make conclusions about the value of H. In particular, in this case 
we have 

V^'^ « b^T + V2b^Tn-^''^N(Q, 1), 

so the behaviour of V^'^ is essentially the same as that of the quadratic power 
variation of Wiener process 

n-l 
k=0 

This confirms that we cannot make any conclusions about H from the asymp- 
totic behaviour of quadratic variation. 

Nevertheless, we will study the behaviour of quadratic variation in this case 
in more detail in order to be able to distinguish between the cases H < 3/4 and 
H > 3/4: statistically. 

Define 

Proposition 3.8. For H G (3/4,1), the sequence {Zf^^ Z}^^\^ • ■ ■ ) converges 
in distribution as k oo to a sequence (CijC2;---) of independent standard 
Gaussian variables. 

Remark 3.4. We emphasize a sharp contrast with the case H e (1/2,3/4), 
where the sequence {Zk,k > 1} has a positive limit in view of (26), hence, 
it eventually becomes positive. This clearly gives a possibility to distinguish 
statistically between cases H e (1/2,3/4) and H e (3/4,1). (See 4.1.3 for 
comparative simulations.) 

Proof. Define 

0^/2 n — k/2 2^ — 1 / 2 \ 

&^^(V'i-T)^^g(2'(Ar..-) -t)^ 

By CLT, ^ ^(0, 1), k oo, so we need to study the collective behaviour. 
To this end, observe that the vector (^^,^^+1, . . . ,^fc+m) can be represented as 
a sum of independent vectors 

2''-! 

{£,k, £,k+l, ■ ■ ■ , £,k+m) — ^ Ck,i, 
i=0 

where the jth coordinate of (^k,i, j ~ 0,1,2, ... ,m, is 
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(We simply group terms on the intervals of the partition 1^2^*^, i = 0, . . . , 2*^}.) 
Therefore, we can apply a vector CLT and deduce that for every m > the vec- 
tor (^fe, ^fe+i, . . . , S,k+m) converges in distribution to an (m-|- l)-dimensional cen- 
tered Gaussian vector as A; — >■ oo. Consequently, the sequence (^^,^^+1,^^+2, . . . ) 
converges to a centered stationary Gaussian sequence as fc — > 00. 

Taking into account the fact that the asymptotic behaviour of V^^'^ and 
is the same (up to the factor &^), we obtain that 

9fc/2 

with U^'^ is defined in (25), also converges to a centered stationary Gaussian 
sequence. It is straightforward to check that the limit covariance is that of an 
i.i.d. standard Gaussian sequence, whence the result follows. □ 



3.14. H = 2,/A 

For H — 3/4, an analogue of Proposition 3.8 can be proved, that is, (Zj., Z^^i, . . . ) 
converges in distribution as fc — > cx) to a sequence (<^i, C2, ■ • ■ ) of independent 
Gaussian variables with unit variance. However, since 

y^3/4,2 ^ j^2j. ^ ^2^3/2^-1/2 ^ ^^2 j,^-l/2^(o^ 

the limiting stationary distribution now has a positive mean, namely, E [ ^1 ] = 
a^&^^T^/^(l — 2^^/^). As long as this value depends on how big is a compared 
to 6, we might be unable to distinguish this case from H > 3/4. On the other 
hand, if b is small relative to a, it might be hard to distinguish this case from 
H < 3/4. 



3.2. Statistical estimation using 4th power variation 

It was mentioned in the previous section that the performance of quadratic 
variation estimators in the case H G (1/2,3/4) is not very satisfactory. One 
could try to improve it by considering quartic variation of 



fc=0 i=0 

Theorem 3.1. Statistics 



Tjii) ^ 1 ttHA 



and 

Hk -7. l0g2+ — ffj 

are strongly consistent estimators of the Hurst parameter H G (1/2, 3/4) in the 
mixed model (15). 
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Remark 3.5. Both these estimators arc quite poor. Regression of several values 
of log2 U^'^ on k leads to a much better estimator. However, as numerical 
experiments in Section 4 suggest, it is better to use the quadratic variation 
based estimators (although they are not apparently very efficient as well). 

Remark 3.6. It is clear that (28) also allows to construct approximate confidence 
intervals for H. One should be careful with the third term in (28): for H G 
(5/8,3/4) it is smaller than the Gaussian fluctuation, so can be omitted, while 
for H E (1/2, 5/8] it does affect a confidence interval. However, a priori a value 
of H is not known, of course, so in case where the value of estimator is close to 
5/8 one should take a partition sufficiently fine in order to determine whether 
if > 5/8 or not. 

Proof. By LLN and CLT, 

V^'^'° w 3T^n-^ + AV6T^n-^/^N{0, 1). 
Using (7), we obtain 

yH,3A ^ T"+'/^^JI^n-"-^N{0, 1) = \/l5r^+3/2„-i-^iV(0, 1). 
Further, from (4) 

yH,2.2 ^ j.2H+l^-2H ^ ^^2^^^^2 + & +1^-2H -1/2 ^ ^ 

Finally, 

Kf -i'^ « ^r3^+i/2n-3«7V(o, 1) 

by (7) and 

-O'^ « 3T4»ni-4H + T^Hn'/^-^"aHAN{0, 1), 

by the ergodic theorem and (11). Collecting all the terms and omitting the terms 
of lower order, we get 

V,^'^ « ST^^n-^ + (ST^^+^a'h'n-^" + 3T*"n^-^"a^ + iVeTH^n-'/^NiO, 1). 

(28) 

Hence, the assertion follows. □ 

3. 3. Estimate of Hurst parameter for known a and b 

When the scale coefficients a and h are known, the estimation procedure signif- 
icantly simplifies, and the quality of estimators is improved. Since we already 
know asymptotic expansions of power variations, we immediately proceed to 
statistical conclusions. 
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3.3.1. He (0,1/2) 

Proposition 3.9. For H e (0, 1/2) and known a, the statistic 

2_ fc + 2 log2 a - log2 V"'"^ 

2(fc-log2 T) 

is a strongly consistent estimator of H . 

Using the asymptotic expansion (24), we can construct another estimator 
of H , which has an asymptotically better rate of convergence, namely, we can 
define Hk{a) as a solution from the interval (0, 1/2) of the equation 

log2(-C/f'2) = 21og2a + 2i7fe(a)log2T + log2 (l - 2i-2ff.(a)^ + ^(^ _ 2i?,(a)). 

We can also construct asymptotic confidence interval using the knowledge of 
a and b; there is no big advantage of this knowledge in the case H S (0, 1/4), 
but for H € (1/4, 1/2) the asymptotic confidence interval will be better. 

Proposition 3.10. Let He {0,1/2), a, b are known. The interval [i/min(a, &), ^^max(a, b)], 
where 

i^min(«, b) = ^{l+ l0g2+ Yk(b) - l0g2+ (>fc+l(5) + Z^<j'^^^^^V^f^2-'^/^)) , 

i/,„ax(a, b) = H^Ua, V ^ (l + log2+ Yk{b) - log2+ (Vk+iib) ~ z^a'^^^^^V^,^^2-'''^)) , 

Yk{b) = V^,'^ -b'T. 
is an asymptotic confidence interval for H of level \ — a. 
Proof. We have 

2^-^HyH,2 _ yHa^ ^ ^^i-2H „ ^ a^T^" cj'h2''^^/^-^"^ N [Q, 1), 

therefore by Slutsky's theorem, 

2'^/'- ^-^ ^^^^^4^ — ^ N{0, 1), fc ^ oo, (29) 

as required. □ 

3.3.2. H e (1/2,3/4) 

In this case the knowledge of a and b allows us we can improve on estimation 
of H by utilising (17). 

Proposition 3.11. If H E (1/2,3/4) and a, b are known, the statistic 
^ k + 2\og^ a -log J V^,'^ -b^T 



is a strongly consistent estimator of H . 
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If only a is known, we can still use (26) to construct an estimator for H . This 
requires solving a non-linear equation 

log2 C/f = 2 log2 a + 2Hf ) (a) \og^ T + (1 - 2i?f ^ {a))k + log^ (l - 2^-^"^' 

on the interval (1/2,3/4). 

Knowing coefficient b also allows to construct more efficient confidence inter- 
vals (those constructed using (27) are pretty useless). 

Proposition 3.12. Let H ^ (1/2,3/4:), b is known. The interval [Ha^inib), Haiii^{b)], 
where 

HmUh) = ^ (l + log2+ Mb) - log2+ (Bu{b) + ^Bk(b)^-Ak{b)Ck{b))) , 

HmUb) - i?mi„(6) V ^ (l + log2+ Ak{b) - log2+ (Sfc(fe) - VSfe(6)2- Afe(6)Cfc(5))) , 
Afe(6) = rfe(6)2 - 4Dk{b), Bk{b) = 2Yk{b)Yk+iib) - Dk{b), 
Ck{b) = Yk+i{b)^ - 2Dk{b), Dk{b) = zlb^T-'2-^-\ 
is an asymptotic confidence interval for H of level 1 — a. 
Proof. Write 

2^--^"(vf,^'-b'T)^(vf,:,\-b'T) 

5TV23-4ff - 2i-2ff 1 ^ ^ ^ ^ 

Then the required result follows by solving the quadratic inequality for 2^^^^ 
in 

□ 

4. Simulations 

In each procedure we take T — a = b = l, n~ 2"^^ and use the circulant 
method to simulate values of on the uniform partition {iT/n, i = 0,1, . . . ,n} 
of [0, T]. The simulation takes about 6 seconds on Intel Core i5-3210M processor, 
computing all estimators takes 10-500 milliseconds. 

4.I. Estimators based on quadratic variation 

4.1.1. He (0,1/2) 

The following table contains estimators of Hurst parameter H and 99 % ap- 
proximate confidence interval [i?min , ^^max] for it, the values of H range from 
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0.05 to 0.5 with step 0.05. For each value of H we also estimate coefBcient a, 
and for H between 0.25 and 0.5 we give value of estimator for b. {N/A means 
that 6|q is negative; * means theoretical value.) 



H 


0.05 


0.1 


0.15 


0.2 


0.25 


0.3 


0.35 


0.4 


0.45 




.046 


.092 


.1382 


.1841 


.2302 


.276 


.3215 


.3656 


.4054 


H20 


.0509 


.1017 


.1489 


.1993 


.2511 


.3018 


.3525 


.407 


.4605 




.0566 


.106 


.1449 


.1922 


.2547 


.2919 


.3534 


.375 


.4397 




.0448 


.0932 


.1311 


.1768 


.2371 


.2708 


.3269 


.335 


.3551 


-^^max 


.0686 


.119 


.159 


.208 


.2728 


.3136 


.3809 


.417 


.5* 


020 


1.019 


1.023 


.985 


.991 


1.014 


1.027 


1.043 


1.135 


1.293 


ho 












2.627 


N/A 


1.823 


1.241 



We see from the table that the estimator H). has consistently the best perfor- 
mance: the relative error is within 2 % for all values but 0.45, where the relative 
error is 2.33 %. 

The estimator underestimates all values of H by around 8 %, which seems 

~ (2) 

a lot for such number of computations. Finally, the relative error of H^. ranges 

from 1 % for H = 0.35 to 6.25 % for H = 0.4 and 13.2 % for H = 0.05. 

— (2) 

This hints at a large variance of the estimator ' , which is supported by the 
following figure. It contains values of estimator Hk (lozenges), Hk (squares) and 
- 5,...,20, = 0.3. The values of Sj^^ outside [0.1,0.5] 



^k"' (triangles), k 



(2) 

were dropped. It is clearly visible that the performance of Hj, is better than 
that of Hk, despite the former has asymptotically better rate of convergence to 
H. 



0.5 
0.45 

0.4 
0.35 

0.3 



0.2^ 
t 

0.2 
0.15 
0.1 





A 



A 





A 






□ g 



□ □ 

A 





□ □ □ 



10 



15 



A 

_0_ 



6 □ o ^ 



20 
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For sake of completeness, we present some values of estimators for H > 0.5. 



H 


0.5 


0.6 


0.7 


0.8 


0.9 


H20 


.4353 


.4577 


.4601 


.4603 


.4604 


H20 


.5021 


.507 


.5027 


.4996 


.4997 



As expected, the values of H20 are close to 0.5. 
4.1.2. He (1/2,3/4) 

The following tabic contains estimators -^20 ^^"^^ -^20 '^^ Hurst parameter H, 
and estimator 620 of the coefficient b. The values of H range from 0.525 to 0.725 
with step 0.025. We also give a "regression" estimator H^^\ It is obtained in the 
following way: we consider linear regression of {log2^ U^'^ ,j = m, m+1, . . . , 20} 
on {to, to + 1, ... , 20}, where m — 11, 12, . . . , 15, and take the best regression 
(in terms of the coefficient of determination). If f^^^ is the coefficient of the 
best linear regression, we set = (1 _ f(2))/2. We do not give values of 

the estimator 0,20 and confidence intervals in view of their poor quality. (In the 
table, N/A means an estimator value outside the theoretical range (1/2,3/4).) 



H 


0.525 


0.55 


0.575 


0.6 


0.625 


0.65 


0.675 


0.7 


0.725 


&(2) 
-"20 


.6061 


.604 


.6096 


.627 


.639 


.6447 


.6629 


.7013 


.6931 


5(2) 

-"20 


N/A 


.5559 


.5341 


.6308 


N/A 


N/A 


.7241 


N/A 


N/A 


5(2) 


.5668 


.5797 


.6083 


.6041 


.6116 


N/A 


N/A 


N/A 


N/A 


&20 


1.229 


1.123 


1.065 


1.033 


1.017 


1.009 


1.005 


1.002 


1.001 



As we can see, when the value of estimator H^^' is greater than .625, it is quite 
reliable, while other two estimators are not. In other cases H^'^^ performs slightly 
better. 



4.1.3. He (1/2,3/4) 

Below are values of {[10'^U"'%k = 11, 12, ... , 20} for H = 0.7 and H = 0.8. 
The difference is clearly visible: for H = 0.7 the sequence is positive, while for 
H = 0.8 there is a plenty of negative values. 



= 0.7 


869 


649 


523 


3 


260 


18 


78 


98 


53 


50 


if = 0.8 


665 


-620 


482 


-475 


8 


-29 


-104 


-71 


-78 


-28 
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4-2. Estimators based on quartic variation 

4.2.1. He (1/2,3/4) 

The following table contains estimators and of Hurst parameter ff, 
the values of H range from 0.525 to 0.725 with step 0.025. We also give a 
"regression" estimator i/^*). It is obtained in the following way: we consider 
linear regression of {\og2+ Uf''^,j = m,m + 1, . . . , 20} on {m, m + 1, . . . , 20}, 
where m = 11, 12, . . . , 17, and take the best regression (in terms of the coefficient 
of determination). If f^"*) is the coefficient of the best linear regression, we set 
if = — f^**) /2. {N/A means an estimator value outside the theoretical range 
(1/2,3/4).) 



H 


0.525 


0.55 


0.575 


0.6 


0.625 


0.65 


0.675 


0.7 


0.725 


-"20 


N/A 


N/A 


.5103 


.5224 


.5303 


.5644 


N/A 


N/A 


.5802 




N/A 


.6227 


.6098 


.5736 


N/A 


N/A 


N/A 


N/A 


N/A 




.5321 


.5606 


.6189 


.6783 


.6378 


.7425 


N/A 


N/A 


N/A 



We see that the estimators based on the quartic variation are quite useless and 
definitely worse than those based on the quadratic variation; only the regression 
estimator is reliable for H close to 0.5. 



4-3. Estimation when a and b are known 

4.3.1. He (0,1/2) 

Below are values of estimators H20 and H20 for H from 0.05 to 0.45 with the step 
0.05. Since the estimates are very close to the real values, we give a rescaled 
difference. We can see that i?20 h^is more stable performance, while H20 has 
worse (but still very good) quality for H close to 1/2. 



H 


0.05 


0.1 


0.15 


0.2 


0.25 


0.3 


0.35 


0.4 


0.45 


10HH2oia)-H) 


.66 


-2.5 


2.3 


-1.3 


.36 


19 


190 


300 


960 


10HH2oia)-H) 


-220 


-22 


51 


50 


-70 


60 


11 


1.1 


47 



Further we give asymptotic 99 % confidence intervals for H e [1/4, 1/2) using 
(29). 



H 


0.25 


0.3 


0.35 


0.4 


0.45 


Hminia,b) 


.2469 


.2972 


.3473 


.3956 


.4457 




.2538 


.3037 


.3537 


.402 


.453 



We see that the intervals are smaller than those given by (3.4). 

The picture below illustrates the quality of estimators Hk{a) (lozenges) and 
Ffe(a) (squares), /c = 5, . . . , 20, iJ = 0.3. 
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0.32 



0.31 



0.3 



0.29 



0.28 



O 
□ 



O 



-Q — ^ — a 



10 



15 



20 



4.3.2. H G (1/2,3/4) 

The following table contains estimators {a, h) and {a) of Hurst param- 
eter H . The values of H range from 0.525 to 0.725 with step 0.025. We also 
give 99 % confidence intervals constructed with the help of (30). Starting from 
H = 0.65, these confidence intervals already contain whole values from 1/2 to 
3/4, so we give 80 % confidence intervals (marked by ^). As before, by * we 
denote the theoretical bound. 



H 


0.525 


0.55 


0.575 


0.6 


0.625 


0.65 


0.675 


0.7 


0.725 


Hil\a,b) 


.5249 


.55004 


.57496 


.5995 


.6262 


.6487 


.6736 


.6927 


.712 




.5538 


.5547 


.5734 


.5962 


.6336 


.6422 


.6642 


.6701 


.683 




.5203 


.539 


.5586 


.571 


.5404 


.609t 


.5951^ 


.6251^ 


.5785^ 


-ffmax(&) 


.5301 


.557 


.5944 


.639 


.6785 


.732t 


.75* 


.75* 


.75* 



. — ("21 — (2) 

We see that {a,b) outperforms (a), but the advantage of the latter is 
that it uses only knowledge of a. The confidence intervals are acceptable for 
iJ < 0.6 and pretty useless for bigger values. 
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